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We review the exact solutions in modified gravity. It is one of the main problems of mathematical 
physics for the gravity theory. One can obtain an exact solution if the field equations reduce to 
a system of ordinary differential equations. In this paper we consider a number of exact solutions 
obtained by the method of separation of variables. Some applications to Cosmology and BH entropy 
are briefly mentioned. 



I. INTRODUCTION 



- - Gravitational field equations describing the geometry of space-time play a fundamental role in modern tlieoretical 

- physics. Their analysis is an extremely difficult task. However, one can find an exact solution in some cases, imposing 
certain additional restrictions. There are several ways to impose constraints on the space-time, for example, the 
, algebraic classification of the Wcyl tensor (of Pctrov types) and the Ricci tensor (Plcbanski type), the choice of 
energy-momentum tensor from physical considerations, presence of the symmetry groups acting on a manifold, etc. 
^ , One can get the exact solution by reducing the system of equations to a system of ordinary differential equations. 
This can be done using the method of separation of variables. In essence, the separation of variables is the only 
I currently known method of structural integration of the field equations. The purpose of the method consists of 
classification of all the privileged coordinate systems and external fields, which is a separation of variables. In the 
Q classification we refer to the transfer of all relevant space-time metrics (non-equivalent with respect to admissible 
coordinate transformations) satisfying the requirement of complete separation of variables in the equations of motion 
of test particles. In flat space-time classification is carried out completely. 
(^JQI From a mathematical point of view, the study of homogeneous spaces and Stackel metrics is of special interest. 
The spaces are united by existence in the space of sets consisting of three geometric objects (for Stackel spaces — 
Killing vector and tensor fields, for homogeneous spaces — Killing vector fields). For both cases the field equations 
' ^ ' '^^^ ^® reduced to a system of finite (but sufficiently large) number of ordinary differential equations. There are 
' many methods in mathematical physics for studying such systems of equations. For example, one can use methods of 
. additional symmetries of the system of equations, the Hamiltonian formulation of the theory of dynamical systems, etc. 
■ Essentially, all physically interesting cases (FRW cosmology, BH) belong to Stackel and homogeneous spaces. 

p: 

O " II- EXACT SOLUTIONS IN STACKE AND HOMOGENEOUS SPACES 

(N" 

Recall that metric is called the Stackel one if the Hamilton- Jacobi equation 

5"^5„5^j=m2 = (1) 



X: 
^3: 



can be integrated by the method of complete separation of variables. The privileged co-ordinate set exists for 
which complete integral of Equation ([1} can be shown in the form 

n 

s = Y,M^\x) (2) 

i=l 

where is the essential parameter. 

It appears that the other important equations of motion (Klcin-Gordon-Fock, Dirac, Weyl) can be integrated by 
complete separation of variables only for the metrics belonging to the class of Stackel spaces. 

That is why the research of this class of spaces belongs to the one of the important problems of the mathematical 
physics. 

The metrics of the Stackel spaces can be used for integrating the field equations of General Relativity (GR) and 
other theories of Gravity. Note that such famous solutions as Schwarzschild, Kerr, NUT, Friedman and others belong 
to the class of Stackel spaces. Apparently the first papers devoted to the problem of classification of the Stackel spaces 
satisfying the Einstein equations were published by Carter [l| . Later in our paper Q the complete classification of the 
special Stackel electrovacuum spaces has been found. In other words, all Stackel spaces satisfying the Einstein-Maxwell 
equations for the case when potentials Ai admit complete separation of variables for Hamilton- Jacobi equation have 
been found. In our paper Q the classification problem has been solved for the case when Ai arc arbitrary functions 
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and spaces are null (types N.l). In our paper (j] all electrovacuum spacetimes admitting diagonalization and complete 
separation of variables for the Dirac-Fock-Ivanenko equation were found. 

One of the complicated problems of the modern mathematical physics is the integration problem of the Einstein- 
Dirac equations. 

Using spaces for which Equation ([l} can be integrated by the complete separation of variables and separated 
solutions of the Dirac equation, one can transform Einstein-Dirac equations to the set of functional equations. The 
first papers devoted to the classification problem for the Einstcin-Dirac equations were done by Bagrov, Obukhov, 
Sakhapov Q. The Stackel spaces of type (3.1) for Einstcin-Dirac and Einstein- Weyl equations have been studied. 
Appropriate solutions have been obtained. They contain arbitrary functions depending on null variable only. 

The problem of classification of Stackel spaces for other theories of gravity for the first time was considered in 
papers 

We have solved the classification problem for the Einstein- Vaidya equations. Let the stress-energy tensor have the 
form 



Tai3 = a{x)lali3, IJ" = 



(3) 



The solution of this problems, as well as a detailed overview of the theory of Stackel spaces, can be found in [lO|. 

Of high interest are homogeneous spaces, which lie at the heart of modern cosmology [lH-[i3|- The homogeneous 
spaces are a base for building the Big Bang model, the initial singularity, as well as the inflationary model. It is of 
interest to identify the various mechanisms of isotropization of the universe [T5j-[20|. Homogeneous spaces are also 
used in a variety of modern theories of gravity for the study of general regularities in the picture of the universe [21 1- 
[2^ . One can study effects of the gravitational field to other fields and matter on the background of homogeneous 
spaces m-H. 

Let us consider the Einstein-Weyl equations. One can show that for all types of Bianchi classification, a closed 
self-consistent system of ordinary integrable diff'erential equations can be constructed [33]-[3l|. For example, consider 
the first type of classification of Bianchi [s^. This is one of the most simple cases; to construct a general solution for 
all types is not possible. 

All calculations will be carried out in the Newman-Penrose formalism. Einstein equation takes the form 



' $00 


= GTqoo'o' 


$01 


= GTooo'i' 


$02 


= GTooi'i' 


< $11 


~ GToio'i' + 


$12 


= GToii'i' 


$22 


= GTiii'i' 


A 


1 G 

= 6^-12^^ 



AB' 
AB' 



(4) 



Here H is cosmological constant, G is corresponds to the gravitational constant, $Qb are Ricci spinors and Taba'B' 
is energy-momentum tensor, 



Tab'cd' — 'ik{^D>V AB'^C + S.b'^cd'O 



^C^AB'£,D' — £,A^CD'£,B' ~ VD'^ AB'VC 
Vb'^CD'VA + Vc'^AB'TlD' + TjA^CD'TlB') 



(5) 



and Vab is spinor derivative. 
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Ricci spinors are expressed in terms of spin factors as follows 

$00 = Dp — Sk — p^ — aa — p(e + e) + Ttr + k(3q! + /3 — tt), 
$01 = -DS — Se — al]} + e — 2e) — /3<7 + pe + kX + kj — tt (e + p) , 
$02 = -DA — Slf + A(3e — e — p) — Jla — tt (tt — /3 + X) + Vk, 

$11 = i [D7 - Ae + (5q; - - Q;(r + tt) - ^(t + tt) + 7(e + e) 

+ s{"f + 7) — TTT + Ki^ — /.jp + aX — aa — fiji + 2a/3— 

= 7(p-p)-e(M-M)], 
$12 = (57 — A/3 — 7(r — a — /3) — pr + crj/ + eI7 + /3(7 — 7 — yu) — aA, 
$22 = (Ji' — A/.i — p^ — AA — /i(7 + 7) + Ftt — v(t — 3/3 — a), 
6A = 2[Jt - Ap - pjl- (7X + t(I3 - a-r) + ^(7 + 7 + i^k]- 

= 1)7 + Ae + - ^ + a{T + 7f) + ;3(r + tt) - 7(e + e)- 

= £(7 + ^) + TTT — KU — pp + X<7 

Choose a metric in the form 

500 = 1, goa = 0, gij = -7jj 

where 7^- is the metric of a three-dimensional space with the signature (+,+,+). It is simple to establish that this 
space permits a three-parameter Abelian group of motions, and hence is a Steckel space of type (3.0). 
The orthogonal tetrad is chosen in the form 

e(o)a = (1,0,0,0), e(i)„ = (0,A,B,C) 
e(2)a = (0, K, S, V), e(3)„ ^ (0, P, M, Z) 

where A, B, C, R, S, V, P, M, Z are arbitrary functions of the time. 
Using this tetrad, we construct the Newman-Penrose tetrad 

mi = --j=(e(^2)i + «e(3)j), = —{e(^2)x - ie{3)i) 
and obtain the following relations between the spin factors 

A = — CT, V ~ K, TT ~T, 7 = —e 

a = /3, p = -p = JI, a = i (r - k) 

The spinor field of spatial rotation of the tetrad may be diagonalized and made real (one real component remains). 
For this case, the energy- momentum tensor takes the form ([SJ. 

f2/-= 



Tqo'oo' = 2zA:<o(e - e), Too'oi' = ika^^r, T^vav = ~2ika^^Qa 

^Ol'lO' 



loi'ii' = ikaS,QK, Tii'ii' = 0, Too'ii' = ika^l{e - e), T( 



The field equations take the form 



(6) 
(7) 
(8) 

2e~ p + p^ + aa -2pe - Akk + 4e^ = — (9) 

(10) 
(11) 



a a — 


2pe- 


4kk 


= 


a — 


2pa^ 


-2ae 


= 2iy^o^cr 


k — 


2pK ^ 


- 2k£ 


= *yCo^K 


2pe - 




h4e2 


H 

~ Y 


P - 


-IpN 


- 2p£ 


H 




He- 


p)io 


= 



Multiplying Equation (O by cr and Equation ([8]) by k, we add and subtract the resulting equations and their 
conjugate forms. Transformation yields 



Finally, we obtain the following equation 



^2 3if 



(1) Consider the case when H = 0. In this case 

ie-p) 



9 = -2 



and thus 

1 



Incorporating the constant into the definition of x, we write 



~ 2x + const 



2x 

Hence 



2x 2 u) X 
where d is a constant of integration. We now find P and e in explicit form 



X X 



2x 2x X X 

where / is a constant of integration. The constants of integration satisfy the condition 

Obviously, this equation has a nontrivial solution. 

(2) Similarly, for the case H ^ 0, we can show that it is not implemented. 
To find the tetrad wc have to solve the two equations 



Since 



(12) 



(InKK) = — 4(e — p) 
{Inaa) = — 4(e — p) 

Multiplying Equation pip by j/^Qj we obtain 

[In(y^o')] =-2(e-p) (13) 

One find that 



(14) 



{e-p) + 2ie-pr = — (15) 



(e-p)^^ (16) 



-- LJ = - (17) 



p+-=0^p=- (18) 



1 1 / / + 1/2 



2Ae - (K - iP)K - {K + iP)K = A' (20) 
2Ak + {K - iP)a + (X + iP)p = -{K' + iP') (21) 

^^1/2 + /^ p=/, K^k-e^-, ^ = 5-e'(2.+c) 

XXX X 



they take the form 



XXX 

2A—e''^ + (K - iP)—e'^^^+^\K + iP)^-^ - (K' + iP') 

XX X 

We multiply the second equation by e~*'^ and make the substitution 

X = K cos u! + P sin uj 
Y = P cos a; — sin w 



The inverse transformations are 



K = Y cos UJ + X sin uj 
P = X cos ui — Y sin oj 



We also introduce the notation ^ = F, ^ = 

A{2F + D)- k{X - iY) - k{X + iY) = 

2Ak + s{X - iY)e"' + F{X + iY) = 
We make one more substitution 



A^ 
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- {X + iY)' -iuj'iX + iY) 



A = We-^'^, {X + iY) = {U + iV)e- 



and take into account that 



A' dA dx dA 
uj' dx duj duj 



Then the equations take the form 

W{3F + D)-k[{U + iV) + iU- iV)] 



dW 
duj 



2kW + se'^iU -iV)^ --^iU + iV) - i{U + iV) 

dw 



If e'" ^ a + ip, I = 3F + D and + (3^ = 1, then 



Wl - 2kU = 



dW 

duj 



2kW + C2aU + {sP -l)V = - ^ 

djUj 

{sp + l)U - saV = - ^ 
duj 



One can obtain a third-order equation in V 

P^-P^l- ?^(«' + 4fc' - 1) + V{4kha + l{s' - 1)) = 
aw^ dw^ duj 

The characteristic equation takes the form 

- X^l - X{s^ + 4fc2 - 1) + (Ak^sa + l{s^ - 1)) = 

From this it is possible to find three roots Ai, A2, A3 

A3 - X^l - A(s2 + 4fc2 - 1) + {APsa + l{s'^ - 1)) = 
A^ - A2(Ai + A2 + A3) - A(AiA2 + A1A3 + A2A3) - A1A2A3 = (4fc2sa + l{s^ 

Ai + A2 + A3 = / 
A1A2 + A1A3 + A2A3 = - (s^ + 4fc2 - 1) 

A1A2A3 = - (4fc2sa + /(s2- 1)) 
In all, four cases are possible: 
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1. Ai, X2, A3 are pairwise unequal and real. Then 

V = cie^i" 



C2e ^ + cae ^ 



2. Ai is real, A2 = A3 



F = cie^i" + C2e^^'^+C3e^^'^ 



3. Ai, A2 = A3 arc real 



V = cie^i" + (C2 + csuj)e 



4. Ai = A2 = A3 = A 



V = (ci +C2a; + C3a;2)e^'^ 
U 



W 



2a(6/3 + 1) 



Thus, an accurate solution of the Einstein- Weyl equation has been obtained for space-time of Bianchi type 1. The 
solution only exists when H ^ 0. The spin coefhcients and the Weyl spinor take the form 



X 



1/2 + / 



sd 



= — e , y^o = - 

X X 



where /, fc, s, y, c are constants and a; = dlnx 

On integrating the equations for the tetrad, the functions obtained will depend on the form of Equation (P^ . As 
an example, the functions for the first case, when all the roots are different and real, are 



A = e- 



cie^i"(Ai" - P + 1) + C2e^^'"(A2' - &^ + 1) + C3e^^"(A3" - P + 1) 

2a(6/3 + l) 

x,u, {ba-M) A,^ (&a-A2) A.^ l^^-Ag) 
' 6/3 + 1 ^ ' 6/3+1 ^ ' 6/? + l 



where Ai, A2, A3 are the roots of Equation (|22p : _F, a, 6, ci, C2, C3, a, /3 are constant. 

The other functions specifying the tetrad are of analogous form. The only difference is that another set of arbitrary 
constants must be chosen in place of ci, C2, C3. 



III. MODIFIED GRAVITY 



In the previous sections a few exact solutions of the classical theory of gravitation were obtained. However, in recent 
years, modified gravity theory is very popular. These theories my be studied also using the methods described above. 
Typically, solutions with a spatially fiat metric depend only on time. It corresponds to the first type of Bianchi and 
type (3.1) on the classification of Stackel. 



The most popular models are models of modified gravity |33j - |35| . which represents a classical generalization 
of general relativity (modifications of the Hilbert-Einstein action by introducing different functions of the Ricci 
scalar [SJ] - [s^ or Gauss-Bonnet invariant H^-ji^), should consistently describe the early-time infiation and late- 
time acceleration, without the introduction of any other dark component. 

In the framework of these theories, a number of cosmological models have been constructed to adequately describe 
the current observational data [4^-[5l|. In addition, under this theory there may be objects such as black holes. The 
properties of these objects are different from classical ones. For example the black hole entropy in the model F (R) 
gravity will have the form (HH, [H^l 

S=^f\R,) (23) 



where Ah = Airr^^. 

As an ex amp le of the exact solution in modified theory of gravity. We consider the sixth dimensional Gauss-Bonnet 
theory [13, [HH • 

We shall start from the following string-inspired action in six dimensions 



S = j (fx^iR + eLGs) 



where e is a constant and Lqb is Gauss-Bonnet invariant 

Lgb = R^.c.pB!"'''^ - ^Rt.uR^" + (24) 

The metric is the product of the usual metric corresponding to the 4-dinicnsional FRW universe and a 2-dimensional 
surface, namely 

ds' = -dt^ + a\t) [{dx^f + {dxy + [dx^f] + b\t) [{dx^f + {dx'^f] , (25) 

the scalar curvature is 

„ 6a2 12d6 262 4^ 
ab 0'' a 

while the four-dimensional and topologically invariant Gauss-Bonnet Lagrangian, Lgb, has the form 

_ iSa^b 72d^b^ 2\b?a 96dd6 24d62 ^%a% 48d66 

Lgs — — 37— ^ 212" ^ 3 ' 27~ uT ^ tT ^ UT^ 

a'^b a'^b'^ w'b ab'' a'^b ab'^ 

or, equivalently, 

24 . . . . ..... 

Lgb = -^i2a^bb + iaa^b^ + a^'db^ + Aaddbb + a^db^ + 2ad^bb + 2a^dbb). (28) 

The corresponding equations of motion are obtained by variation of the action with respect to a and 6, which yields 

d^b^ + 4ad56 + a^b^ + 2adb^ + 20^})}) + \2ed^b^ + 16edd56 4- Seadb^ + Sed'^bb + IGeadbb = 0, 

3ad'^b + 3a'^db + 3a'^db + a^b+l2ed^b+12ed'^db + 2Aeaddb+l2eaa^b =0. (29) 

These equations can be easily rewritten in terms of the Hubble rates H — a' /a and h — b' jb^ namely 



3/1^ + 4/iF + 3i72 + 2h + 2H + l&eh^H + 28eh^H^ + IGehH^ + IGehhH + 8ehH^ + 8eh^H + IGehHH = 0, 
+ 3hH + 6H^ + h + 3H + 12eh^H^ + 36ehH^ + 12eH^ + 12ehH^ + 2AehHH + UeH^H = 0. 



(30) 



In addition, variation over the metric in the above expressions gives the constraint equation 

+ 6hH + 3ij2 + 36eh^H^ + 2AehH^ = 0. (31) 
This equation helps to exclude h and h' from Equation pop . As a result, one gets an equation for H only: 

H' -3H^ X ^ + "^"^ + <-22^/6 + 6AG)H'^ - 24e2(9V6 - 52G)H^ + 96e3(17V6 + 12G)H^ - SOGAV&e'^H^ 
^ V6-12e(V6 + 16G)i?2 + 72e2(3\/6-32G)i/4- 2880V6e3i?6_^ 31104 Vee^ifS 

where 



G = \J\ - 6em + 24e2iJ4 (33) 
One can check that this last equation obeys the fundamental relation (for e > 0): 

H (^2_^2)(^2_,2)/W (34) 

where p, q and r are constants, and the function f{H) < 0. 
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For e = one does recover (as it should be) an explicit solution. In other cases, only a numerical analysis can be 
carried out. In fact, the equations for a, b are 

b^a'^ + 4aba'b' + aH'^ + 2ab^a" + 2a%b" = 0, iaba'"^ + ia^a'b' + ia%a" + a%" = (35) 

and, from here, 

3/1^ + AhH + 3i/2 + 2h' + 2H' = 

(36) 

+ 3hH + 6H^ + h' + 3H' = 

from where one gets that i/' = (3 ± 2^/6)H^ and the solution is given by if = — ^jj^^ being a = 3 ± 2\/6- Moreover, 
in terms of the scale factors: 

a = C2[±(at + Ci)]-i/" 
6 = C3[±(at + Ci)]^/" 

where /3 = 3±-\/6. General approach to Stackcl spaces of first sections may be applied to such theory as well. However, 
this is more complicated in modified gravity. 



IV. CONCLUSIONS 



In this article we consider the exact solutions constructed in the classical theory of Stackel spaces and partially 
homogeneous spaces. We consider in detail the exact solution of Einstein- Weyl space of the first type of classification 
of Bianchi in the Newman-Penrose formalism. As an example we consider the exact solution for the modified gravity 
theory of the Gauss-Bonnet for the six-dimensional metric which depends only on time and has a diagonal form. 
Thus, we discuss the problem of obtaining exact solutions in the different theories of gravity. The paper shows how, 
by use of the method of separation of variables, can one construct exact solutions of cosmological models, both for 
space with matter (the classical theory of gravity) and in the case of the modified theories of gravity (six-dimensional 
theory of the Gauss-Bonnet). 

Note that separation of field equations for modified gravity in Stackel space may be done in analogy with the method 
developed at the beginning of this work. However, the corresponding generalization is very cumbersome technically. 
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